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We argue that the critical dynamical fluctuations predicted by the mode-coupling theory (mct) 
of glasses provide a natural mechanism to explain the breakdown of the Stokes-Einstein relation. 
' This breakdown, observed numerically and experimentally in a region where MCT should hold, is 

one of the major difficulty of the theory, for which we propose a natural resolution based on the 
recent interpretation of the MCT transition as a bona fide critical point with a diverging length scale. 
^ We also show that the upper critical dimension of MCT is dc = 8. 

o 

^ Mode Coupling Theory (mct) provides a useful theoretical framework to account for many of the empirical proper- 

ties of liquids, at least in the weakly supercooled region 0,0 ■ There are however a number of well known difficulties 
with the theory, the most prominent being the absence, in real systems, of the dynamical arrest singularity predicted 
by MCT. This very singularity underlies most of the quantitative predictions of MCT, such as the divergence of the 
, relaxation time, the power-law shape of the correlator in the [3 region, or the critical behaviour of the non-ergodic 
(Edwards- Anderson) parameter q Additional activated relaxation processes, not described by MCT, have to be 
invoked to argue that the MCT transition temperature should be understood as a cross-over rather than a true 
G singularity. In the region between Tc and the so-called onset temperature, Tonset, where the slow dynamics regime 
sets in 3j, MCT fares quite well to describe a host of experimental and numerical results with one noticeable excep- 
r-| , tion: the breakdown of the Stokes-Einstein relation. This breakdown is one of the most important aspect of the 
Q ' phenomenology of supercooled liquid. The Stokes-Einstein relation states that the product of the viscosity 77 times 
O , the self-diffusion constant D, divided by the temperature is a constant independent of temperature. This relation is 
well obeyed in liquids at high enough temperature. Instead, in supercooled liquids Drj/T increases quite shaply below 
I \ Tonset :4| and eventually reaches values of the order of 1000 or more at the glass temperature Tg In other words, 
^ the self-diffusion of particles becomes much faster than structural relaxation. The decoupling between viscosity and 
10 diffusion is interpreted by many to be a direct piece of evidence for dynamical heterogeneities in glassy liquids j^. 
. Intuitively, this decoupling comes about because diffusion is dominated by the fastest particles whereas structural 
relaxation is dominated by the slowest regions [5|. Detailed explanations based on different theoretical approaches 
have been put forward in the literature iQiQiHTIS- In most models, the breakdown of the Stokes-Einstein relation 
is in fact a direct measure of the width of the distribution of local relaxation times (see below and H, IE S IS 
, Standard MCT does not account for the sharp increase of Drj/T between Tonset and To although this is precisely the 
"Y^ region where MCT is supposed hold. The aim of this short contribution is to discuss an explicit mechanism that leads 

' to such a breakdown within the mode-coupling theory of glasses, once critical fluctuations are taken into account. 
^ , In its usual interpretation, MCT describes homogeneous dynamics and neglects all fluctuations; it therefore seems 
I ' that any effect associated to dynamical heterogeneities is outside of the scope of MCT. This pessimistic view was 
y however recently argued to be erroneous 

I1SIIII1III3 (see Ullil for earlier insights). MCT should in fact be seen as 
g a standard mean-field theory, except that it provides a self consistent equation for the two-body dynamical correlation 
^ function (describing density fluctuations) rather than for a one body order parameter, such as the magnetisation in 
• • ■ the usual Curie- Weiss theory of magnets. In the latter case, we know that the appearance of a non trivial solution of 
the mean field equation below a certain critical temperature is in fact associated with the divergence of the magnetic 
susceptibility, itself related to magnetisation fluctuations. These fluctuations are harmless in high enough dimensions, 
$H ' but become dominant in dimension space less than four, where all critical properties of the phase transition are 
5^ , strongly affected by these fluctuations. The same scenario, although more involved, also holds for MCT above a 
certain dimension dc, MCT is expected to be quantitatively accurate (at least in a regime where the above mentioned 
activated events can be neglected), whereas for d < dc, dynamical fluctuations play a major role and must be properly 
accounted for. In this case, fluctuations of the (two body) dynamical correlation are measured through a four body 
correlation G4(r, t); its integral over space define, in analogy with ferromagnets, a dynamical susceptibility called 
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X4{t) in the recent literature 0,0,0- The extension of MCT to compute Gi{r,t) reveals that Xi indeed diverges 
as the mode-coupling temperature is approached 0, 0| . This indicates that the dynamics becomes correlated 
over larger and larger length scales as the system freezes, which is in fact expected on general grounds: a diverging 
relaxation time should necessarily involve an infinite number of particles Il9l. The growth of has recently been 
detected in numerical simulations and experimentally in supercooled liquids [T^ llTl l2Cl|| and in granular media |2l| ; 
it should also transpire in the divergence of the non-linear susceptibility of glassy systems . 

As discussed in 'll] the spatial correlations underlying the mode-coupling singularity necessarily lead to the existence 
of an upper critical dimension for this problem. In [Tlj it was shown that in the early /^-relaxation regime, the 
four-point correlator reads, in Fourier space, Gidk, t ~ Tp) ^ (fc^ + y^)^^, where e — {T ~ T^jT^,. From this, one can 
estimate, in the spirit of a Ginzburg argument 11], the fluctuations of the non-ergodic parameter g in a region of size 
^ ~ to find ^'^iq ~ ^ Imposing that Sq must be much smaller than the critical behaviour predicted by MCT, 

i.e. qc — q^ \/e, one finds that this assumption is only consistent when d > dc = 6. For d < 6, the fluctuations become 
dominant for e small enough and change all the critical exponents. But most interestingly, as we will argue below, 
these fluctuations are also responsible for the breakdown of the Stokes-Einstein relation. Before doing so, we should 
however point out that the analysis of 11] in fact overlooked a contribution due to the coupling between dynamic 
fluctuations and slow conserved degrees of freedom such as density and energy fluctuations j2J]. These fluctuations, 
close enough to the transition, become the leading ones and change the above mean-field critical behaviour of G4 
to G4{k,t w Tp) ~ (fc^ -I- \ft)~'^ 23]. In this case, the fluctuations of the non-ergodic parameter in a region of size 
^ grow as ^"^+4/2^ which becomes dominant below the critical dimension dc = 8(38|]. A diagrammatic derivation of 
this upper critical dimension is presented in the Appendix. Note that in reality the situation is more complicated: 
from numerical simulations of Lennard-Jones systems [2,'i ] one finds that fluctuations due to slow conserved degrees 
of freedom only become dominant very close to the transition, where MCT breaks down. For colloids, where MCT 
fares rather well, we do not have yet quantitative estimates of the relative contribution of these fluctuations. It might 
therefore well be that in cases of experimental or numerical interest, the original analysis of BB holds and the effective 
upper critical dimension is dc = 6. In any case, the aim of this paper is not to obtain quantitative predictions but 
only to showing that critical dynamical MCT fluctuations provide a natural scenario for the decoupling of self-diffusion 
and viscosity. 

In principle, the contribution of critical fluctuations to the viscosity and diffusion constant should be calculated 
using a renormalisation group around dc- Here, we do not attempt to do this but argue physically that such a program 
should lead to a decoupling between these two quantities. We assume that non trivial critical exponents describe, in 
d < dc, the different physical quantities; for example ^ ~ e"", qc — q^€^,T^ e and G4{r) ~ g{r /^) /r'^~^^'^ . Mode- 
coupling theory describes how the random potential created by density fluctuations slows down the particles. Close 
to the MCT transition, this self-generated random potential persists on a timescale comparable to the one responsible 
for structural relaxation. Hence, it acquires a static component that self-consistently traps the particles at the MCT 
transition, where the relaxation timescale diverges. The strength of this self-generated component, and hence of the 
trapping potential, is measured by the plateau value q of the correlation function, which therefore acts as an effective 
coupling constant. When q reaches q^, the random potential is sufficiently strong to prevent the particles from moving. 
Clearly, activated effects are expected to destroy the transition (any finite potential barrier can be overcome in finite 
time), but within MCT, these processes are neglected. 

The local difference between q and q^ therefore plays a major role: if the local value of q is slightly smaller, particles 
are less trapped and diffusion is enhanced. If, on the other hand, the local value of q is slightly larger, particles 
are frozengSEj]. In other words, the self-generated disorder fluctuates in space and leads to fluctuations of the critical 
temperature on correlated regions of linear size Since the (local) relaxation time is set by the distance from the 

(local) critical temperature, the dynamics strongly fluctuates from one correlated region to another. This is precisely 
why the fluctuations of the local correlation function, measured by Gi{r,t), play a crucial role to generate strong 
dynamical heterogeneities. Parallel to the Ginzburg argument, the fluctuations (per particle) of the non-ergodic 
parameter in a region of size ^ ~ are given by 5q ~ ^(4-d-i?)/2 ^ ^-v{i-d-ri)/2 ^ compared with qc — q ^ e^ ■ 

These fluctuations of the effective coupling constant in turn induce local fluctuations of the relaxation time as: 

5t Sq _ ^ -(^+^.-4) _^ 

T f3qc-q 

For d > dc, the mean fleld value of the exponents v — 1/4, /? = 1/2 can be used, and one finds that St/t 
close to the MCT transition: dynamical heterogeneities are mild and do not jeopardize the Stokes-Einstein relation, as 
predicted by MCT calculations that neglect spatial fiuctuations altogether. Interestingly, this conclusion does not hold 
for d < dc, where the relative width of the local relaxation time distribution must increase as e ^ 0. For e not too 
small, one can still use the mean-field value of the exponents, but now the fluctuations of local relaxation times are 
found to increase as St/t ~ e-l-ic-rf)/*^ before one enters the Ginzburg region, where fluctuations eventually change 
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the values of the exponents. For standard phase transitions, the hyperscahng relation v{d+r] — 2) — 2(3 = is such that 
fluctuations remain of order one when e ^ 0. Its generalization to the mct transition would be v{d + 77 — 4) — 2/3 = 0. 
However, it is not clear whether within mct hyperscahng holds since the MCT transition is purely dynamical, and 
has a mixed first order/second order nature [qc is non zero ai T — Tc while the relaxation time is diverging as a 
power law). A detailed renormalisation group study of the replica field theory with a cubic term could help settling 
this interesting technical point. In any case, we expect that in three dimensions, neglecting the role of the activated 
processes (that cut off the transition) and entering in the regime where the dynamics is supposed to slow down because 
of the MCT critical point (i.e. between Tonset and Tc) 5t /t should first grow as e decreases, and either saturate or 
diverge as e — > 0, depending on the validity of hyperscahng. 

As already discussed, a key ingredient that generates violation of the Stokes-Einstein relation is the fluctuation of 
the local relation timescale r. Different models capture this effect. For instance, assuming as in local fluctuations 
of the diffusion coefficient one finds that the Stokes-Einstein violation is related to the ratio between (1/r) and 1/(t). 
Another model consists in assuming that particles hop with random relaxation times. In this case the self diffusion 
coefficient is given by -D = -^^/(r), where £ is the typical hopping distance, expected to be of the order of the particle 
size and roughly temperature independent. The viscosity, or the terminal relaxation time, on the other hand, are given 
by the integral of the average correlation function, which weighs slow regions proportionally to the local relaxation 
time; therefore rj ^ (r^) /(r). In both models the crucial ingredient to evaluate violation of the Stokes-Einstein relation 
are indeed local fiuctuations of the relaxation times. For example for the latter model, one gets: 

cxl + ((,5r)2)/(r)2^1 + Ce'^('^+"-4)-2/3^ (2) 



KseT 



where C is a numerical constant of order unity and Kse is the high temperature value of Drj/T. For d < dc, this 
formula predicts an increasing decoupling between viscosity and diffusion between Tonset and Tc. If C is not too 
small, we predict a power-law viscosity-diffusion decoupling as Tc is approached, with exponents that change as one 
enters the Ginzburg region. Such a behaviour may have been observed in numerical simulations of Lennard- Jones 
particles 0, where Dr]/T is found to increase by a factor ~ 5 in the MCT region P|, while MCT predicts an increase 
of a mere 1.2 nil. These results are interesting because they show that both the viscosity and the inverse diffusion 
constant can be made to scale relatively well with e = (T — Tc)/Tc with the same Tc for both quantities, but with 
different exponents. In other words, this scenario is compatible with a 'fractional' Stokes-Einstein relation, D ^ -q^^ 
with fc < 1. More precisely, the numerical results are 77 ^ e"^ '', whereas D ^ e^'^, so that Drj ~ g-O-6 (^ggg Q a,nd refs 
therein) . These results suggest that the observed breakdown of the Stokes-Einstein relation might indeed be related 
to critical fluctuations. A similar breakdown of the Stokes-Einstein relation in the mode coupling region has also 
been observed in hard sphere systems, both experimentally |26j and numerically ^29l ISfljl . Of course, much stronger 
violations are observed closer to Tg, but outside the region where mct can be valid. These should be related to the 
broadness of activated relaxation time distributions and are outside the scope of the present discussion [4lj . 

Summarizing this short contribution, we have argued that a dramatic consequence, in low space dimensions, of 
the critical dynamical fluctuations predicted by the mode-coupling theory of glasses is the breakdown of the Stokes- 
Einstein relation. This breakdown, observed numerically and experimentally in a region where MCT should hold, 
i.e. above Tc, is one of the major difficulty of mct (setting aside other well known difficulties related to activated 
processes), and for which we provide a natural interpretation. We believe that this matter deserves more investigation, 
both 



theoretically, to estimate the value of the exponents describing this breakdown in d < dc and analyze the 
interplay between activated processes and critical dynamical fluctuations; 

and numerically and/or experimentally to ascertain the connection between dynamical fluctuations and viscosity- 
diffusion decoupling in the T > Tc region, for example by comparing the location r and the width St of the 
peak of the four-point susceptibility Xiit). 

Numerically one could also test the dimension dependence of Stokes-Einstein violation, and compare with other 
scenarii, such as suggested by Kinetically Constrained Models which also predict such violations [31|. If conflrmed, 
those effects would be a direct proof of the importance of critical fluctuations in glasses. This would put the program 
of extending mct using renormalisation group methods on top of the agenda. 
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FIG. 1: Diagrammatic representation of the parquet diagrams, k denotes the wave- vector entering into the ladders. 
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FIG. 2: Diagrammatic representation of the squared-parquet diagrams, k denotes the wave-vector entering into the ladders. 



Appendix 

In the following we shall show that below eight dimensions there are infrared divergences that change the mean-field 
scaling derived in BB and |2^. 

The natural framework to analyze critical dynamical fluctuations and their role in determining the upper critical 
dimension is dynamical field theory. The starting point to analyze the dynamics of super-cooled liquids consists in 
writing down some exact or phenomenological stochastic equations for the evolution of the slow conserved degrees 
of freedom. Field-theories are obtained through the Martin-Siggia-Rose-deDominicis-Janssen method, where one 
first introduces response fields enforcing the correct time evolution and then averages over the stochastic noise [ssf . 
Standard mct corresponds to a self-consistent one-loop approximation (see Refs. 0, |3^ |3^ for a discussion of the 
different field-theories and subtleties appearing in the derivation of mct). 

Here we do not provide any detail and we refer to [23Ll33| for a thorough presentation. We only recall that generically 
the four-point function can be written in terms of the so called parquet diagrams, see (371] for a general introduction 
and for an application to supercooled liquids. As shown in (23l | the four-point function can be exactly expressed 
as a sum of parquet and "squared-parquet" diagrams (we will use the same notation of "2^), see Fig. QJand Fig. |21 

The elementary block used to construct the parquet dia gram s corresponds to [5Yj / 6G]GG where S is the self-energy 
considered as a function of the dressed propagator G, see [23 ■ mct consists in retaining only the bubble diagram for 
the self-energy. In this case the elementary block is very simple and leads to the ladder diagrams studied in BB and 
|23l |. which lead to the mean-field scaling discussed in the main text. The diagram in Fig. Ogives a contribution to 
G4(fc,t) that, for times t in the /3 regime and for small k, has a critical behavior like l/(fc^ + -\/e). The diagram in 
Fig. 121 on the other hand, gives a contribution that has a critical behavior like 1/(A:^ -I- y/ef ■ 

From the point of view of critical phenomena, one expects that adding diagrams other than the bubbles to the 
self-energy should be harmless and should not change the scaling for all dimensions larger than the upper critical 
dimension dc- In order to show that this is the case and determine dc one has to prove that (1) adding any type of 
diagram to the bubbles does not change the mct (mean-field) critical properties in presence of an infrared cut-off 
on momenta integration and (2) equating the infrared cut-off to zero leads to a change of the mean-field scaling 
only below the u ppe r critical dimension dc- The first part of this procedure is lengthly and difficult and it will be 
shown elsewhere [S^l- Concerning the second part, we shall only show an example of diagrams that are responsible 
for changing the mean-field scaling below eight dimension. (A full analysis showing that above eight dimensions all 
diagrams are harmless is outside the scope of this article and is left for future work). 

In order to show an example of divergent diagram below eight dimension one has to focus on the elementary block 
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FIG. 3: Divergent diagram below eight dimension assuming mean-field scaling. One box represents the parquet diagram shown 
in Fig. Note that the lines on the left have to be considered without propagators. 

[JS /5G\ GG used to construct the parquets. The small k behavior of this block determines the mct mean- field scaling 
discussed above {k is the momentum entering into the block as shown in Fig. 1 and 2). Within mct, and more 
generally, if one puts an infrared cutoff, one finds a scaling fc^ -I- ^/e at small k. One can show that this assumption is 
inconsistent in low enough dimension, focusing on an explicit diagram, see Fig. ^ |42l|. 

In order to do that, we shall assume that the mean- field scaling is correct and compute the small k behavior of 
this diagram. In dimensions less than dc, the corresponding contribution to [5Ti/5G]GG at small k will be found to 
differ from mean-field. For the diagram in Fig. |31 the small k dependence is given by (we omit the dependence on 
the non-critical wave- vectors): 



For e = the small k behavior is then given by /(O, 0)-f c/c^-f ... where /(O, 0) and c are two well-defined constants above 
eight dimensions (that can be explicitely computed from /(fc, e)). /(O, 0) and c lead, respectively, to a renormalization of 
the critical temperature and of the "rigidity" constant in front of the fc^ term in G^. However, below eight dimensions, 
c diverges (if the infrared cutoff vanishes) and the behavior of /(fc, 0) is singular at small k: /(O, 0) + c'fc^^'^*^'^^ + .... 
This clearly shows that the assumption of mean field scaling is not correct. In order to go further and compute the non 
mean field exponents one should perform a renormalization group treatment, or some self-consistent (mode-coupling!) 
resummation of these dangerous diagrams. We leave this for further investigations. 

As a conclusion we find, in agreement with the simple Ginzburg argument developed in the main text, that below 
eight dimension critical fluctuations are not governed by mean-field theory. 

Finally, we remark that for a system with no conserved degrees of freedom, as it is the case for the Langevin 
dynamics of disordered p-spin systems, the squared-parquet diagrams are absent. Therefore the divergent diagram 
corresponding to the one in Fig. Ois now formed by single parquets both on the top and bottom parts of the diagram. 
Repeating the same analysis performed above, one finds dc = 6 in this found in BB. 
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